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and Turbulent Heating 
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An experimental program has been conducted to obtain data on the effects of distributed 
surface roughness on blunt bodies. Sphere-cone and hemisphere geometry wind tunnel 
models with both pattern and sand-grain distributed surface roughness were tested at Mach 
6 across a free stream Reynolds number range of Reoo = 3.0xl0 6 /ft to Re^ = 8.3xl0 6 /ft. 
Laminar, transitional and turbulent surface heating data were obtained using global 
phosphor thermography. Roughness-augmented, turbulent heating levels from 1.0 to 1.9 
times the predicted smooth-wall turbulent levels were measured. 

Nomenclature 

= average roughness element spacing (in.) 

= pattern roughness element diameter (in.) 

= heat-transfer film coefficient (kg/m/s 2 ) 

= heat-transfer film coefficient based on Fay-Riddell theory (kg/m/s 2 ) 

= heat-transfer film coefficient based on smooth wall, turbulent prediction (kg/m/s 2 ) 

= heat-transfer film coefficient based on experimental roughness data (kg/m/s 2 ) 

= tunnel total enthalpy (J/kg) 

= wall enthalpy (J/kg) 

= enthalpy at 300 K temperature (J/kg) 

= physical roughness height (in.) 

= equivalent sand-grain roughness height (in.) 

= Reynolds number based on friction velocity and roughness height 
= pattern roughness element spacing (in.) 

= free stream Mach number 
= free stream pressure (Pa) 

= heat-transfer rate (W/cm 2 ) 

= heat-transfer rate based on Fay-Riddell theory (W/cm 2 ) 

= radial coordinate (m) 

= base (max) radius of model (in. or m) 

= corner (shoulder) radius of model (in. or m) 

= nose radius of model (in. or m) 

= free stream unit Reynolds number (1/m or 1/ft) 

= momentum thickness Reynolds number 
= roughness Reynolds number 
= free stream temperature (K) 

= boundary-layer edge velocity (m/s) 

= velocity at roughness height (m/s) 

= free stream velocity (m/s) 

= friction velocity (m/s) 

= Cartesian coordinates (m) 

= angle of attack (deg) 
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= boundary-layer height (m) 

= sphere-cone model cone angle (deg) 

= momentum thickness (m) 

= Dirling spacing parameter (in.) 

= boundary-layer edge density (kg/m 3 ) 

= density at roughness height (kg/m 3 ) 

= density at wall (kg/m 3 ) 

= free stream density (kg/m 3 ) 

= boundary-layer edge viscosity (kg/m/s) 
= wall viscosity (kg/m/s) 

= shear stress (Pa) 


I. Background and Introduction 

In this work, the effects of surface roughness on boundary-layer transition and turbulent heating are considered. 
These effects are important because the surface roughness of an entry vehicle’s Thermal Protection System (TPS) 
can promote earlier boundary-layer transition and produce higher turbulent heating (and shear) levels than would be 
expected based on an idealized, smooth- surface analysis. However, due to the complexities of roughness effects, a 
vehicle’s TPS is typically designed using analytical, computational and/or experimental techniques that are based on 
the assumption of an Outer Mold Line (OML) with a smooth surface. The effects of roughness on the 
aerothermodynamic environment are then included through approximate engineering correlations and methods. 

“Roughness” is a generic term in aerospace literature that can encompass different types of surface features, as 
illustrated in Table 1. Roughness can be divided into two general types: discrete and distributed. Discrete 
roughness includes surface features such as attachment-point cavities, compression pads, gaps or steps between 
heat-shield tiles and physical damage to a TPS. Distributed roughness includes features such as tightly- spaced TPS 
tile patterns or hexcombs and ablated TPS. 

The historical database on discrete and distributed roughness effects (e.g. surveys presented in Refs. 1 and 2) has 
largely been built on the basis of ground testing of simple geometries such as hemispheres, flat-plates, and sphere- 
cones and from flight data on small-angle cones (e.g. missiles) and lifting bodies (e.g. the Space Shuttle). There are 
less data (especially for distributed roughness) available on blunt-body entry vehicles such as typified by the Apollo 
(spherical cap segment), Mars Science Laboratory (70-deg sphere-cone) or Genesis (60-deg sphere-cone) 
configurations. 

To address the lack of distributed roughness data on blunt bodies, an experimental program has been conducted 
to obtain hypersonic wind tunnel data on the effects of distributed roughness on blunt bodies. Two blunt-body 
geometries were considered in the study: a 70-deg sphere-cone and a hemisphere. Two types of distributed 
roughness were employed: sand-grain and patterned. Data were obtained on these geometries to determine the 
effects of roughness on boundary-layer transition and turbulent heating augmentation and are presented herein. 


Table 1. Surface roughness types 
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II. Test Program 


A. Experimental Facility 

Testing was performed in the NASA Langley Research Center (LaRC) 20-Inch Mach 6 Air Tunnel (Figure 1). 
This tunnel is a blow-down facility in which heated, dried, and filtered air is used as the test gas. A detailed 
description of this facility can be found in Ref. 3. The tunnel has a two-dimensional contoured nozzle that opens 
into a 20.5 in. x 20.0 in. (0.52 m x 0.508 m) test section. The tunnel is equipped with a bottom-mounted injection 
system that can transfer a model from the sheltered model box to the tunnel centerline in less than 0.5 sec. Run 
times of up to 15 minutes are possible in this facility, although for aeroheating studies run times of only a few 
seconds are typically required. The nominal reservoir conditions of this facility produce perfect-gas free-stream 
flows with Mach numbers between 5.8 and 6.1 and unit Reynolds numbers of Reoo= 0.5><10 6 /ft to 8.3><10 6 /ft 
(1.64xl0 6 /mto 27.4xl0 6 /m). 
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Figure 1. LaRC 20-Inch Mach 6 Air Tunnel 


B. Data Acquisition, Analysis and Uncertainty 

1. Phosphor Thermography 

Convective aeroheating measurements were made using the two-color, relative-intensity, global thermographic 
phosphor method (Refs. 4-5). In this method, a ceramic model coated with a thermographic phosphor compound is 
illuminated by ultra-violet light sources that produce temperature-dependent fluorescence of the phosphor coating. 
Images of the model are taken in the tunnel before and during a run using a three-color, charge-coupled device 
camera. The fluorescent intensity data are converted to temperatures based on camera-response calibrations. Heat- 
transfer rates are then determined with the assumption of a one-dimensional, step-function heating response (which 
corresponds to a parabolic temperature-time history) beginning at injection of the model into the tunnel. 

The IHEAT (Imaging for Hypersonic Experimental Aero thermodynamic Testing) code is used to perform the 
data reduction. IHEAT provides global heat-transfer distributions expressed in terms of the ratio Mz FR , where /z FR is 
the heat-transfer coefficient resulting from a Fay-Riddell computation (Ref. 6) for the heat-transfer to a sphere of 
user- specified radius at cold- wall (300 K) conditions. Data are typically output as 2-D image files. Quantitative 
values can also be extracted from the image data along specific reference lines with geometric (x-y) parameters that 
are approximately corrected through linear interpolation between user- specified reference locations. 

2. Three-Dimensional Data Mapping 

The image data produced by IHEAT are corrected for optical perspective effects and mapped to a three- 
dimensional (3-D) surface model of the test article. To accomplish this mapping, iterative perspective, rotational, 
and translational transformations are performed on the 3-D surface geometry until its 2-D projection matches that of 
the 2-D image data. The image data are then assigned transformed (jc, y, z) coordinates based on interpolation 
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between the image and surface geometry, and then the transformation is inverted to obtain a 3-D heating 
distribution. After mapping the image data to the 3-D surface geometry, direct comparison can be made to 
computational heating predictions. 

3. Data Uncertainty 

For testing of blunt-body configurations such as the in this study, an uncertainty estimate for the phosphor 
thermography measurement technique of approximately ±10% can be made based on Ref. 5. Uncertainties also are 
introduced through the variation of wind-tunnel flow conditions. Since the repeatability of flow conditions is very 
good, the resulting uncertainty in heating measurements is estimated to be only on the order of ±5%. The other 
major source of uncertainty that must be considered is the mapping of 2-D image data to a 3-D surface. This 
uncertainty is highly problem dependent; on the overall acreage of a blunt-body geometry, this uncertainty is 
estimated to be on the order of ±10%. However, in geometrically-complex regions, such as around small-scale 
features where the scale of such a feature approaches the resolution of the camera optics, this uncertainty could be 
higher. For this work, the acreage uncertainty of 10% is applied and resulting total uncertainty based on a root- 
mean-square summation of acreage mapping, flow conditions and measurement technique gives an estimate of 
±15%. 

C. Roughness Model Design and Fabrication 

Three types of wind tunnel models were tested in this study: smooth- surface models, models with patterned 
surface roughness elements, and models with sand-grain surface roughness. Smooth models, pattern roughness and 
sand-grain roughness models were fabricated for the sphere-cone geometry, and smooth and sand-grain roughness 
models were fabricated for the hemisphere geometry. The pattern and sand-grain roughness sphere-cone models are 
shown in Figure 2 and Figure 3 and the hemisphere sand-grain roughness models are shown in Figure 4. The 
sphere-cone model diameter was 8 -in. and the hemisphere diameter was 6-in; other geometry parameters are listed 
in Table 2. The model designation and fabrication process for each type of model is discussed below. 


Table 2. Model geometry parameters 


Model 

R B 

R B 

R N 

R N 

Rn/ Rb 

(in) R( (m) 

R C f /?b 

6c 


(in.) 

(m) 

(in) 

(m) 




(deg) 

Sphere-cone 

4.00 

0.1016 

2.00 

0.0508 

0.500 

0.200 0.00508 

.005 

70 

Hemisphere 

3.00 

0.0762 

3.00 

0.0762 

1.000 

NA NA 

NA 

NA 


i 10x15 I 10x30l | 20x301 i 230-Mesh^ Bl 40-Mesh IT 80-Meshl 

CuulBlLj 



Figure 2. Pattern roughness sphere-cone models Figure 3. Sand-grain roughness sphere-cone models 
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Figure 4. Sand-grain roughness hemisphere models 


1. Smooth-surface wind tunnel models 

Wind tunnel models for use with the phosphor thermography method are slip-cast in silica-ceramic and coated 
with a thermographic phosphor compound. The standard fabrication process (Ref. 7) begins with the production of 
a rapid-prototype, resin pattern of the geometry. The resin pattern is built in a stereo-lithographic apparatus (SLA) 
from a CAD data file. These patterns are then hand-worked to a smooth surface finish to remove any SLA 
manufacturing artifacts. A multiple-piece injection mold (from which the resin model can easily be removed) is 
then built around the resin pattern. The resin pattern is then removed and wax is injected into the mold to form a 
wax pattern. Next, a new two-piece shell mold is built around the wax pattern, and then the wax is burned out of the 
shell mold. A silica ceramic model is then slip-cast in the shell mold. Then, the ceramic model is removed from the 
shell mold, dried and sintered. The finished ceramic model is then back-filled with a hydraulically setting magnesia 
ceramic for strength and support. Finally, the model is coated with a mixture of phosphors that luminesce under 
ultraviolet lighting. This process was used to create smooth- surface models to serve as a baseline for comparison 
with roughness effects. 

As a reference point for comparison with the roughness models, the surface roughness of the phosphor coating 
has been determined in previous measurements (Ref. 8) to be approximately 0.0005 in., which is below the range at 
which roughness effects on transition or heating would be expected. Turbulent data obtained in previous studies 
produced close comparisons between measured and predicted heating levels, which supports the conclusion that the 
phosphor surface roughness does not affect measured heating levels. The phosphor coating process does however, 
occasionally produce small, discrete surface features (on the order of 0.0015 in.) that can affect transition onset, but 
not cause heating augmentation. Similar effects can also sometimes be observed as a result of damage to the coating 
that occurs in small areas during handling and testing of the models. 

2. Patterned surface roughness wind tunnel models 

For the patterned surface models, metallic molds of the heat-shield surface were first fabricated in place of the 
SLA resin patterns. A chemical etching process was then employed to create hemispherical cavities of various 
diameters and spacings in the molds. Wax patterns were then created from the metallic molds and the rest of the 
fabrication process followed the same steps as for a standard, smooth- surface model. A metallic mold, ceramic 
model, and closeups of the roughness different roughness elements are shown in Figure 5. 

The pattern roughness layout, also illustrated in Figure 5, consists of a repeated hexagonal array of hemispheres 
over the entire heat-shield surface of a model. The defining parameters of a pattern are the hemispherical element 
diameter, D, and the element center-to-center spacing L. Six pattern layouts were considered and the parameters are 
given in Table 3. The pattern name format is “DD x LL”, where “DD” and “LL” are element diameter and spacing 
(both in inches) multiplied by 1000 and the wind tunnel models are identified by this pattern name. 
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Figure 5. Pattern Roughness Layout 


Table 3. Patterned surface roughness model parameters 


Pattern 

Element 

X 

(in.) 

d AV 

(in.) 

D 

(in.) 

L 

(in.) 

k 

(in.) 

^ws 

(in 2 ) 

A P 

(in 2 ) 

k s /k 

k s 

(in.) 

10x30 

hemispherical 

14.07 

0.0279 

0.0100 

0.030 

0.0050 

7.854E-05 

3.927E-05 

0.91 

0.0046 

10 x 15 

hemispherical 

7.03 

0.0140 

0.0100 

0.015 

0.0050 

7.854E-05 

3.927E-05 

3.41 

0.0171 

40 x 120 hemispherical 

14.07 

0.1117 

0.0400 

0.120 

0.0200 

1.257E-03 

6.283E-04 

0.91 

0.0183 

20x30 

hemispherical 

7.03 

0.0279 

0.0200 

0.030 

0.0100 

3.142E-04 

1.571E-04 

3.41 

0.0341 

30x45 

hemispherical 

7.03 

0.0419 

0.0300 

0.045 

0.0150 

7.069E-04 

3.534E-04 

3.41 

0.0512 

40x60 

hemispherical 

7.03 

0.0558 

0.0400 

0.040 

0.0200 

1.257E-03 

6.283E-04 

3.41 

0.0683 


The equivalent sand-grain roughness heights, which are also listed in Table 3, are based on the theory developed 
by Dirling (Ref 9) to correlate pattern-type roughness data with sand-grain roughness data. The equivalent sand- 
grain roughness k s is defined as 

k s 0.0164 x A 378 , for A <4.93 
} k | 139 x A (_1 90) , for A > 4.93 

where the physical roughness k is taken as DU (the particle radius) and the Dirling spacing parameter A is given 
by 



The wetted surface area A ws and profile area A P for a hemispherical element are given by 
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(5) d AV = , 
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The functional relationship between equivalent sand-grain roughness and the Dirling spacing parameter is 
plotted in Figure 6. The resulting values for the sphere-cone models are also shown. 



Figure 6. Equivalent sand-grain roughness 
and Dirling spacing parameters 


3. Sand-grain surface roughness wind tunnel models 

Sand-grain surface roughness was produced by applying an adhesive coating to the SLA resin patterns and then 
dipping the patterns into a bed of precision-manufactured, spherical glass particles. The particles adhered to the 
surface of the patterns to create a distributed, sand-grain roughness over the entire pattern (on the hemispheres, the 
roughness was truncated approximately one inch from the end of the pattern to ensure the roughness elements did 
not cause binding when removing the model from the mold). The patterns were then hand- worked to remove any 
obvious surface irregularities (i.e. clumps of glass spheres) and then ceramic models were fabricated from them 
following the same steps as for smooth- surface models. A resin pattern coated with particles, a finished ceramic 
model, and closeups of the various surface roughness elements are shown in Figure 7. 

The spherical glass particles applied to the models were procured with diameters ranging from 0.0025 in. to 
0.0689 in., where the particle sizes were specified according to an ATSM El 1:01 standard wire mesh standard 
(defined by the diameter of a particle which will fit through a mesh sieve). Models are identified according to the 
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ATSM designation of the particle size and a listing of the models and roughness parameters is given in Table 4. The 
“smooth-wall” baseline model roughness, as defined by the roughness of the thermographic phosphor coating is also 
listed. For the scope of the present study, the sand-grain roughness of a model is simply assumed to be equal to the 
particle radius. While alternative definitions based on statistical analyses of measured surface roughness height 
distributions can be found in the literature, the range of roughness considered in this study spans two orders of 
magnitude and as such it is assumed that the precise definition of the roughness height is a 2 nd -order consideration. 



Figure 7. Sand-grain Roughness Layout 


Table 4. Sand grain surface roughness model param eters 

ATSM Element Diam. Diam. k s 
ID (mm) (in) (in) 

10-Mesh hemispherical 1.7500 0.0689 0.0344 

20-Mesh hemispherical 0.8530 0.0336 0.0168 

40-Mesh hemispherical 0.4250 0.0167 0.0084 

80-Mesh hemispherical 0.1800 0.0071 0.0035 

140-Mesh hemispherical 0.1060 0.0042 0.0021 

230-Mesh hemispherical 0.0630 0.0025 0.0012 

Phosphor hemispherical 0.0254 0.0010 0.0005 


III. Computational Method 


A. LAURA Flow-Field Solver 

Flow field predictions were performed using the LAURA (Langley Aerothermodynamic Upwind Relaxation 
Algorithm) code (Ref. 10) to obtain heat- transfer rates for comparisons with the experimental data and boundary- 
layer parameters for use in correlating the data. LAURA is a three-dimensional, finite-volume solver that includes 
perfect-gas and non-equilibrium chemistry options, a variety of turbulence models, and ablation and radiative 
transport capabilities. 
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In this study, the perfect-gas air model was used for the wind tunnel predictions. Both laminar and turbulent 
solutions were generated, with the Cebeci-Smith algebraic model (Ref. 11) being used for the turbulent cases. 
Solutions were generated on structured, finite- volume multi-block computational grids with a singularity- free nose 
blocks. Grid adaption to the flow features was performed to align the outer boundary of the grid with the shock and 
to cluster grid cells near the surface to produce wall cell Reynolds numbers on the order of 1 . 

Free- stream conditions were set equal to the tunnel operating conditions, with the wall temperature specified as 
a constant 300 K. The use of a constant wall temperature is acceptable because the heat-transfer coefficient varies 
only very slightly over the range of wall temperatures produced in this facility. 

B. Computational Validation 

In this study, values for turbulent roughness heating augmentation are defined based on comparisons of the 
experimental data to computed turbulent heating levels. It is therefore important to establish the accuracy of the 
computational technique. Relevant comparisons to laminar and turbulent smooth- wall sphere-cone heating data 
have been performed previously for studies in the LaRC 20-Inch Mach 6 Air Tunnel (Ref. 12) and the Arnold 
Engineering Development Center (AEDC) Hypervelocity Tunnel 9 (Ref. 13). For laminar conditions, LAURA 
predictions agreed with experimental data to within the estimated experimental uncertainty, with the exception of 
small regions around the stagnation point on the windward flank of the heat shield. For turbulent conditions, 
LAURA predictions were within or slightly below the experimental uncertainty range depending on test conditions. 
This under-prediction of smooth-wall turbulent cases from past studies is important to consider in comparisons to 
roughness-augmented turbulent data from this study. 

A likely explanation for at least part of these differences is in the specification of the transition onset location. In 
LAURA, the algebraic turbulence models do not have an inherent transition onset mechanism. A transition onset 
location can be specified by the user, but lacking a validated and accepted transition criteria for large-cone angle, 
blunt body flows, the entire flow field is usually assumed to be turbulent, as was the case for this study. However, 
as was shown in the LaRC test (Ref. 12), the turbulent heating levels on the flank of a large cone-angle blunt body 
can be dependent on the transition onset location; in general, the farther outboard transition onset is specified to 
occur, the higher the predicted downstream turbulent heating levels. Thus, if a solution is generated assuming fully - 
turbulent flow over the entire vehicle, the heating levels will be lower than measured levels for a case where 
transition occurs downstream of the stagnation point on the conical flank. Presumably this effect is due to the 
boundary layer not having sufficient running length to reach an equilibrated turbulent profile before reaching the 
shoulder. This assumption can be supported by noting that comparisons for the AEDC test (Ref. 13) were better at 
higher free stream Reynolds number - i.e. for cases where the transition onset location was closer to the stagnation 
point. 


IV. Results and Analysis 


A. Test conditions 

Roughness data were obtained in two wind tunnel entries, identified as Test 6971 for the sphere-cone models and 
Test 6975 for the hemisphere models. Runs were performed at free stream Reynolds numbers of Reoo = 3.0><10 6 /ft 
to Re^ = 8.3><10 6 /ft with 16-deg angle-of-attack for the sphere-cone and 0-deg angle-of-attack for the hemisphere. 
Tree stream conditions for the tests are given in Table 5. 


Table 5. Nominal Conditions for LaRC 20-Inch Mach 6 Air Tunnel Tests 6791 and 6975 


Re» 

( 1 /ft) 

Reo, 

( 1 /m) 


Poo 

(Pa) 

Too 

(K) 

Poo 

(kg/m 3 ) 

Uoo 

(m/s) 

H 0 -H 300 K 

(J/kg) 

Hemisphere Sphere-cone 

hpR-hemi QFR-hemi ^FR-hcmi QFR-hemi 

(kg/m 2 /s) (W/cm 2 ) (kg/m 2 /s) (W/cm 2 ) 

3.03E+06 

9.95E+06 

5.99 

844.1 

62.5 

4.708E-02 

948.7 

2.120E+05 

2.163E-01 

4.581 

2.649E-01 

5.611 

5.04E+06 

1.65E+07 

6.02 

1421.0 

63.2 

7.843E-02 

957.5 

2.213E+05 

2.823E-01 

6.239 

3.458E-01 

7.641 

6.63E+06 

2.18E+07 

6.04 

1833.0 

62.6 

1.023E-01 

954.6 

2.170E+05 

3.212E-01 

6.984 

3.934E-01 

8.554 

7.46E+06 

2.45E+07 

6.04 

2058.0 

62.5 

1.154E-01 

953.6 

2.167E+05 

3.402E-01 

7.364 

4.167E-01 

9.019 

8.34E+06 

2.74E+07 

6.03 

2090.8 

58.6 

1.249E-01 

918.1 

1.791E+05 

3.407E-01 

6.101 

4.173E-01 

7.472 
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B. Heating Data and Comparisons with Predictions 

Data were obtained for each sand-grain and pattern roughness element configuration for each model at all free 
stream Reynolds number conditions, plus a small number of repeat runs, for a total of more than 100 separate cases. 
As this data set is too large to be presented in its entirety here, data will be presented for each roughness 
configuration only at the lowest and highest Reynolds numbers of Reoo = 3. Ox 10 6 /ft and 8.3xl0 6 /ft. Analyses and 
comparisons presented are based only on centerline values. As these values represent only a small subset of the 
large, global data set, the results may not be representative of behavior over the entire geometry. 

Global aeroheating distributions from the thermographic phosphor measurements and comparisons of centerline 
data with LAURA laminar and fully-turbulent predictions are shown in Figure 8 and Figure 9 for the sphere-cone 
with sand-grain roughness, Figure 10 and Figure 11 for the sphere-cone with pattern roughness, and Figure 12 and 
Figure 13 for the hemisphere with sand-grain roughness. In all cases, the images and line-plots are ordered in terms 
of increasing sand-grain height, or the for pattern models, by equivalent sand-grain height. 

The observed behavior of the data for both the sphere-cone and hemisphere sand-grain roughness cases is 
consistent and follows expected trends: as the sand-grain roughness height is increased for a given Reynolds 

number, the transition onset location moves closer to the stagnation point and the turbulent heating levels increase 
further over the predicted turbulent values. 

The sphere-cone pattern roughness data behavior is also generally consistent and follows expected trends. 
However, there is one discrepancy in that the heating data from 40x120 pattern model ( k s = 0.0183 in.) shows higher 
augmentation than that of the 20x30 pattern model ( k s = 0.0341 in.), when according to the theory presented for 
equivalent sand-grain roughness heights, the reverse would be expected. With respect to this observation, several 
points should be considered: 1) The Dirling spacing value for the 20x30 model is close to the inflection point in the 
functional relationship between X and k s /k and some uncertainty in the definition there might be expected; 2) other 
authors have suggested alternative forms for this correlation (e.g. Ref. 14) that may provide a better representation 
of these data; and 3) the physical roughness height of the 40x120 pattern is greater than the boundary layer height 
(i.e. k/b > 1) in a significant region between the stagnation point and the nose of the sphere-cone, and thus 
correlations based on the assumption of roughness elements submerged within the boundary layer may not 
accurately reflect the current data set. 
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10-Mesh, Run 81 
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LAURA, turbulent 



Centerline Heating Distribution 


Figure 8. Sphere-cone sand-grain roughness data, Reoo = 3.0 x 10 6 /ft, a = 16-deg 


Smooth, Run 7 230-Mesh, Run 1 


140-Mesh, Run 18 80-Mesh, Run 23 


40-Mesh, Run 28 20-Mesh, Run 33 



0 0.5 1.0 1.5 2.0 






2.0 


Test 6971, Run 7, Smooth OML 
Test 6971, Run 13, 230-Mesh 
Test 6971, Run 18, 140-Mesh 
Test 6971, Run 23, 80-Mesh 
Test 6971, Run 28, 40-Mesh 
Test 6971, Run 33, 20-Mesh 
Test 6971, Run 80, 10-Mesh 



Centerline Heating Distribution 


Figure 9. Sphere-cone sand-grain roughness data, Reoo = 8.3 x 10 6 /ft, a = 16-deg 
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Figure 10. Sphere-cone pattern roughness data, Reoo = 3.0 x 10 6 /ft, a = 16-deg 
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Figure 11. Sphere-cone pattern roughness data, Reoo = 8.3 x 10 6 /ft, a = 16-deg 
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Figure 12. Hemisphere sand-grain roughness data, Reoo = 3.0 x 10 6 /ft, a = 0-deg 
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Figure 13. Hemisphere sand-grain roughness data, Reoo = 8.3 x 10 6 /ft, a = 0-deg 
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C. Turbulent Roughness Heating Augmentation 

A comparison of the current results with data and correlations from previous studies can be made in terms of the 
relationship between the roughness Reynolds number k+ and turbulent roughness heating augmentation, 
trough/ h smooth • The roughness Reynolds number is defined in terms of the friction velocity as: 

(6) k + - P " U ^ - l W T u./A„ 

^ Mw/Pw 

Centerline k + distributions from the LAURA solutions (based on the sand-grain roughness or equivalent) are 
shown in Figure 14 - Figure 15 for the sphere-cone with sand-grain roughness, Figure 16 - Figure 17 for the sphere- 
cone with pattern roughness, and Figure 18 - Figure 19 for the hemisphere with sand-grain roughness. As 
demonstrated by these figures, the test conditions for this study produced roughness Reynolds numbers spanning 
two orders of magnitude, from k + ~ 1 to 100. It has been shown (Ref. 15) for incompressible flat-plat flow - and the 
theory is generally accepted at least qualitatively for compressible flow - that the effects of roughness on the skin- 
friction (and by analogy, heat-transfer) can be divided into three regimes based on the roughness Reynolds number, 
each of which is encompassed by the conditions of the current study: 

1) For k + < ~10, the flow is hydraulically smooth. Roughness elements are fully embedded in the laminar sub- 
layer and have minimal affect on the skin-friction 

2) For ~10 < k + < ~70, the flow is transitionally rough. Roughness elements begin to contribute to skin- 
friction 

3) For ~70 < k + , the flow is fully roughness dominated. The laminar sub-layer disappears, shear stress is 
directly transmitted to the surface by pressure drag on the roughness elements and skin-friction becomes 
independent of viscosity. 

The definition of the roughness heating augmentation is more subjective. Ideally, the denominator in this term 
would be based on smooth-wall experimental data at fully-turbulent conditions on the geometry of interest. 
However, the conditions necessary for smooth-wall fully-turbulent flow are beyond the range of the LaRC 20-Inch 
Mach 6 Air Tunnel; previous experience from testing at AEDC Hypervelocity Tunnel 9 suggests that a Reynolds 
number in excess of Reoo = 30><10 6 /ft would be required. Thus, the predicted fully-turbulent heating levels are used 
in the denominator in place of experimental values. 

^ rough _ ( ^FR ) measured 

^ smooth (h/h 

FR ) turbulent-CFD 

As noted previously, the smooth-wall, fully-turbulent CFD predictions using LAURA with the Cebeci-Smith 
turbulence model are generally lower than measured smooth-wall, turbulent values (especially for cases where the 
flow field is not fully-turbulent over the entire body) which is may be due to lack of a transition model in the solver 
and the assumption of fully-turbulent flow over the entire body. Thus, the turbulent roughness heating augmentation 
definition given by Eq. (7) will introduce a small positive bias to the augmentation values. 

Turbulent roughness heating augmentation values were computed for a single, centerline point on each 
geometry: at z/R = 0.5 (i.e. on the leeward cone flank) for the sphere-cone and at z/R = 0.7 for the hemisphere. 
These locations were selected because the measured heating distributions suggested they were far enough 
downstream from the transition location to be in the turbulent flow regime. These values are presented in Figure 20 
for the sand-grain roughness sphere-cone, in Figure 21 for the pattern roughness sphere-cone, and Figure 22 for the 
sand-grain roughness hemisphere; also, all the data points are shown together in Figure 23. For clarity, data for 
laminar and transition conditions (for which the augmentation ratio was much less than 1, are not shown. 
Additionally, a shaded region is shown on each figure to illustrate an approximate range of historical data sets and 
correlations, such as from Ref. 16. 

For the sphere-cone geometry, heating augmentation levels from approximately h rough /h sm0 oth = 1.0 to 1.8 were 
observed for sand-grain roughness and from h rough /h sm0 oth = 1.0 to 1.6 for pattern roughness. For the hemisphere 
geometry, heating augmentation levels between h rough /h sm0 oth = 1.0 to 1.9 were observed. These data exhibit a first- 
order, logarithmic relationship between roughness Reynolds number and turbulent roughness heating augmentation 
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levels such as seen in prior studies. Additionally, while the data exhibit some scatter, they fall within the historical 
range of data and correlations. 

Because of the limitations noted previously with respect to the assumption in the predictions of fully- turbulent 
flow over the entire and the selection of only a single, centerline data point for comparisons rather than values over 
the whole body, these results should be considered more qualitative than quantitative. In theory, a more accurate 
augmentation ratio could be defined by using the current data set to develop a rough-wall boundary-layer transition 
criteria for use in the LAURA code. A transition criteria would then allow the onset of transition and the 
development of fully-turbulent flow to be included in the flow-field solution, making for a more appropriate 
comparison between data and predictions. And given that the data sets are global, comparisons could be made for 
all points on the body, not just the centerline, to provide a statistically-valid comparison set. However, such work is 
beyond the scope of the initial analysis presented here. 




Figure 14. Sphere-cone sand-grain k+ Figure 15. Sphere-cone sand-grain k+ 

distributions, Reoo = 3.0 x 10 6 /ft distributions, Reoo = 8.3 x 10 6 /ft 




Figure 16. Sphere-cone pattern k+ 
distributions, Reoo = 3.0 x 10 6 /ft 


Figure 17. Sphere-cone pattern k+ 
distributions, Reoo = 8.3 x 10 6 /ft 
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Figure 18. Hemisphere sand-grain k+ 


Figure 19. Hemisphere sand-grain k+ 
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distributions, Re^ = 8.3 x 10 6 /ft 
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Figure 20. Sphere-cone sand-grain roughness 
augmentation (at zJR = 0.5) 
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Figure 21. Sphere-cone pattern roughness 
augmentation (at z/R = 0.5) 
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Figure 22. Hemisphere sand-grain roughness 
augmentation (at zIR = 0.70) 


Figure 23. Roughness augmentation (all data) 


D. Roughness Effects on Transition 

As can be seen in the previously- shown comparisons between predicted and measured heating levels, distributed 
roughness has a direct influence on the location at which boundary-layer transition occurs. The ability to predict 
where distributed roughness causes transition to begin is of value both for vehicle design purposes and for the 
development of turbulent roughness heating augmentation correlations. A correlation for transition onset in terms 
of boundary-layer parameters and roughness height has been given in Ref. 17: 


( 8 ) =- 


500, 

1 0 exp jlog ( 250 ) - log [( Re^ /Re s ) ( ,u w /^ )] } , 
100 , 


for: (Re / /Re fl )(^,/, u ,v) < 0-5 
for: 0.5<(Re t /Re 0 )(/i < ,/,M H ,)<2.5 
for :(Re / /Re fl )(/i ( ,/,«„ ) > 2.5 


where 


(9) Re 0 = ? momentum thickness Reynolds number 

Ve 

and 

(10) Re^ = ? roughness Reynolds number 

AC 

The first region of the correlation corresponds to a “smooth-wall” upper limit for transition and the last region 
corresponds to a rough wall lower limit. Between these two regions, the slope varies linear on a log-log scale and 
this function can be rearranged in terms of the critical transition roughness Reynolds number: 


(11) [Re*], -250 

Centerline transition onset locations were determined at each free stream Reynolds number for all the 
geometries. These values are plotted against the correlation of Eq. (8) in Figure 24 for the sand-grain roughness 
sphere-cone, Figure 25 for the pattern roughness sphere-cone, and Figure 26 for the sand-grain roughness 
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hemisphere; also all data points are plotted together in Figure 27. The data match the general trends of the transition 
correlation, including the upper smooth-wall and lower rough-wall limits and the log-log variation in the roughness 
transition regime. The scatter of the data shown is greater than the uncertainty bounds given in Ref. 17, however 
this initial analysis was limited to only centerline transition. A complete analysis of the transition onset location 
over the entire surface of the models, not just the centerline, will be required to develop a statistically-valid 
correlation for this data set. 




Figure 24. Sphere-cone sand-grain 
roughness transition onset 


Figure 25. Sphere-cone pattern 
roughness transition onset 
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V. Summary 

Roughness-augmented data on boundary-layer transition and turbulent heating have been obtained through 
testing of models with distributed sand-grain and pattern roughness at hypersonic conditions. Testing was 
conducted on sphere-cone and hemisphere geometries at free stream Reynolds numbers of Reoo = 3.0><10 6 /ft to Reoo 
= 8.3><10 6 /ft, which produced laminar, transitional and turbulent flow conditions. The models had physical 
roughness heights from 0.012-in. to 0.034-in, which produced roughness Reynolds numbers spanning two orders of 
magnitude, from k + ~ 1 to 100. This range encompassed flow regimes from roughness having no effects to fully 
roughness-dominated flows. 

Turbulent heating augmentation levels due to roughness exhibited a logarithmic dependency on the roughness 
Reynolds number with values ranging between h rough /h smooth = 1.0 to 1.9, which is consistent with the range of data 
and correlation in the literature. The influence of roughness on transition was analyzed in terms of boundary-layer 
parameters (Ree and Re kw ) and roughness height. The data exhibited similar behavior to classical transition 
correlation, including upper smooth-wall and lower rough-wall limits on transition and a log-log dependency 
between these bounds. 
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